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We present calculations of differential cross sections for Compton scatter- 
ing from the proton, using amplitudes calculated to fourth order in heavy 
baryon chiral perturbation theory. We compare with available data up to 
200 MeV. We find that the agreement for angles below 90° is acceptable over 
the whole energy range, but that at more backward angles the agreement de- 
creases above about 100 MeV, and fails completely above the photoproduction 
^ ■ threshold. 
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I. INTRODUCTION 



Heavy baryon chiral perturbation theory (HBCPT) is a systematic framework in which 
to describe the interactions of nucleons with pions and photons. Respecting gauge and 
Lorentz invariance, it reproduces low-energy theorems, which have been known since the 
50's and 60's, as the first terms in chiral expansions, and in a number of cases also makes 
further parameter-free predictions. One of the simplest and cleanest processes to which it can 
be applied, theoretically at least, is Compton scattering, and indeed this was first studied 
in the very early days of the theory. First spin-independent scattering was considered. 
At lowest order (technically second order in HBCPT) the scattering amplitude is simply 
the energy-independent Thomson term, which is entirely independent of all details of the 
proton's structure. This arises simply gull term in the Lagrangian, one of the terms 

that arise when antinucleons are integrated out of the lowest-order relativistic Lagrangian. 
Energy and angle dependence in the amplitude enter at third order, where the only diagrams 
which contribute are one-pion loop graphs. The next lowest terms in an energy expansion of 
the amplitude (of order u 2 ) have coefficient which are called the polarisabilities, a and (3, and 
at this order in HBCPT they are simply predicted in terms of well-known nucleon and pion 
properties The numerical values, a = 12.5 x 10 -4 fm 3 and (3 = a/10, are in remarkably 
good agreement with the current world averages for the proton, a = 12.1±0.8±0.5x 10 -4 fm 3 
and a + (3 = 14.1 ± 0.5 x 10" 4 fm 3 @. This was an early triumph for HBCPT. 

The calculation of the forward spin-independent scattering amplitude to fourth order 
followed quickly However at that order there is no further predictive power for the 

polarisabilities within the strict framework of the theory, as there are four unknown low- 
energy constants in the fourth-order Lagrangian which give contributions to the isoscalar 
and isovector a and (3. Also, the amplitude is not easily compared with experiment on its 
own. 

The next step was the calculation of the full spin-dependent amplitude to third order p . 
There are four independent amplitudes, each with a leading term (of order u) which is given 
by a low-energy theorem |J. The next terms, of 0(uj 3 ), have coefficients, denoted ji, which 
are called spin polarisabilities. Unfortunately these are not well known experimentally; 
indeed till recently only the contribution which enters forward scattering, called 70, was 
known at all, and that only from pion photoproduction data rather than direct Compton 
scattering experiments [|7|-|9|] . More recently 7^ (for backward scattering) has been extracted 



from unpolarised Compton scattering experiments, with conflicting results [10,11]. The first 



direct measurements of polarised photon-proton scattering at MAMI yield a more reliable 
value 70 of —0.8 ±0.1 x 10 -4 fm 4 Jl2| which is within the range of previous estimates. In 
addition pion photoproduction data has been used as input into fixed-t dispersion relations 
to obtain estimates of the 7$ individually fl~3Hl5l- Unfortunately the lowest order HBCPT 
values are in total disagreement with the experimental numbers, the only saving point of 
agreement being that the isovector quantities, which vanish in LO HBCPT, do seem to be 
smaller in general than the isoscalar ones. 

Most recently, three groups have calculated the spin polarisabilities to fourth order (NLO) 
|i~6Hl9ll- Unlike the spin-independent case, this order contains no LEC's, and so there are 
again genuine predictions. Surprisingly the NLO contributions are as large as the LO ones. 
The agreement with the experimental numbers however is variable. 
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Fig. 1: Diagrams which contribute to spin-dependent forward Compton scattering in the 
e ■ v — gauge at 4th order. The solid dots are vertices from £( 2 \ 

Of course, classic HBCPT has no explicit A. This is assumed to have been integrated out, 
and its legacy is in the LEC's of the nucleon-only theory. It is seen first at different orders in 
different processes: second order for pion scattering, fourth order for spin-independent and 
fifth order for spin dependent Compton scattering. It is expected to be important; it may 
improve the agreement of the spin polarisabilities, but risks spoiling the spin-independent 
ones. 

All this emphasis on the polarisabilities obscures the fact the HBCPT does in fact give 
complete amplitudes for scattering, which one might hope would be believable up to energies 
of the order of the pion mass, since the chiral expansion is both in powers ofu/A and M^/A 
(where A is either m N or the chiral scale 4^/^). When the full third-order amplitudes 
were published, it was claimed that they gave good agreement with cross section data for 
unpolarised scattering up to the photoproduction threshold. However a more wide-ranging 
study was done by Babusci et al. ipOfl , who showed that the agreement started to breakdown 
for backward angles above about 120 MeV. They stopped their comparisons at the threshold. 

We have now calculated the full scattering amplitude to fourth order in HBCPT. Since 
the fourth-order contributions to the spin polarisabilities are so large, it is clearly of interest 
to see whether the full amplitudes and cross sections are also changed dramatically, and 
this is the subject of this paper. We have compared with the available data up to 200 MeV 
21| , |23| -f26| , and find that in fact the difference between third- and fourth-order cross sections 



is small. Above the threshold the good agreement at forward angles continues, while the 
lack of agreement at backward angles is dramatically worse. 
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Fig. 2: The fourth order (solid line), third order (dashes) and Born (dotted) predictions 
for various lab angles. The data are from Baranov, Moscow (green inverted triangles) 



21 1, Ziegler, Mainz (orange stars) |p3[] , Federspiel, Illinois (blue triangles) p4| . |22| , Hallin, 

J, Olmos de Leon, 



Saskatoon (black squares) [p5| , MacGibbon, Saskatoon (red circles) [g 
Mainz (magenta diamonds) |27| and various earlier experiments, listed by Baranov (cyan 
pentagons) [^]. Cross sections are in units of nb/sr, and energies in MeV. Data points have 
been included if their angle is within 2 degrees of the quoted value. 



II. THE CROSS SECTION IN HBCPT 



The usual notation for the scattering amplitude in the Breit frame is, for incoming real 
photons of energy u and momentum q and outgoing real photons of the same energy and 
momentum q', 

T = e^e" 

= e ■ eAi(u;,0) + e ■ qe • q A 2 (u,9) 

+icr ■ (e' x e) A 3 (u, 0) + ia • (q' x q) e' • e A^lu, 9) 
+ (ier ■ (e' x q) e ■ q' — icr ■ (e x q') e' • q^ A 5 (u, 9) 

+ (ier ■ (e x q') e • q' — ia • (e x q) e' • q) A 6 (u, 9), 



(1) 



where hats indicate unit vectors. By crossing symmetry the functions A, are even in uj for 
i = 1, 2 and odd for i = 3 — 6. The amplitudes may be decomposed into pole and non-pole 
pieces. The pole, which occurs at unphysical values of uj and t (= 2uj 2 (cos9 — 1)), arises from 
Born diagrams in which the intermediate nucleon is on shell, and the contribution containing 
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Fig. 3: Predictions for various CM angles. See Fig. 1 for legend. 

this pole can be calculated using Dirac nucleons and on-shell photon-nucleon couplings. The 
leading terms are given by low-energy theorems With the pole terms truncated at fourth 
order (that is, at 0(l/m 3 N )), the amplitudes read 



At{uj,e 

A 3 (u,e 
a 4 (uj,6 

A 5 (lu,6 

Aaiuj.e 



Q 2 e 2 



m f 



Ami, 



((Q + k) 2 (1 + cos6)-Q 2 )(1 — cos 9) uj z 



+An(a + cosd(3) u 2 + 0(uj 4 ) 

k(2Q + k) cosOu 2 - Anf3u 2 + 0(cj 4 



Am 3 N 
e 2 u 
2ml 



Q{Q + 2k) -{Q + k) 2 cos O) + Al + 4vtcj 3 (7 1 + 7s cos 6) + 0(u* 



e 2 oj 



2ml 



{Q + k) 2 + AI +4nu 3 l2 + 0(u* 



e 2 u 



2ml 



(Q + k) 2 +Af+ 47r^ 3 7 4 + 0(u 5 ) 



e 2 oj 



2m 2 N 



Q{Q + k) + Ai + 4vrcu 3 7 3 + 0(uj 5 



(2) 



where the charge of the nucleon is Q = (1 + t 3 )/2 and its anomalous magnetic moment 
is k = (k s + k v t 3 )/2. The polarisability terms are the leading terms in the expansion of 
the full non-pole amplitudes; the polarisabilities are isospin dependent. The contribution 
from the n° graph, A\ , has been separated out, though its contribution is often included in 
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Fig. 4: Effect of different handling of third-order predictions. Solid: prediction taken as 
referring to Breit frame and variable change performed to reproduce the physical value of the 
threshold; dashed: as solid, but no variable change; dotted: predictions taken as referring 
to the centre of mass frame, as in ref. 



the definition of the polarisabilities. Only four of the spin polarisabilities are independent 
since three are related by 75 + 72 + 274 = 0. The spin-dependent amplitudes have no pole 
contributions beyond the leading ones at this order. (At fifth order there are pieces of order 
c<j 3 /m^, though these vanish for both forward and backward scattering amplitudes.) The 
amplitudes are given for a nucleon spinor normalisation of uu = 1. 

To calculate these amplitudes in HBCPT, we need the all the diagrams which appear at 
fourth order. (The pion-loop diagrams are shown in Fig. 1.) There is no clear separation at a 
diagramatic level between pole and non-pole; heavy-baryon seagulls give pieces of the Dirac 
pole but, at fourth order, they also contribute to a and (3. Conversely fourth-order loop 
diagrams have contributions at 0(u) which correct the bare anomalous magnetic moment in 
the LET terms. However, taking all diagrams together HBCPT reproduces the pole terms 
of Eq. with the physical (as opposed to bare) nucleon mass and magnetic moment, as an 
expansion in powers of l/m N . Everything else is non-pole. (While this separation into pole 
and non-pole is useful in discussion, nothing in this paper depends on them. We simply 
calculate the full amplitudes to a given order and use them to produce cross sections. This 
is not affected by the divergence of opinion that is found in the literature over the definition 
of the polarisabilities PSR .) 

For the spin- dependent amplitudes, no LEC's enter except well-known nucleon and pion 
properties. We take g A = 1.267, f w = 92.32, k s = -0.120, k v = 3.706 and M v = 139.6 (as 
the loops are of charged pions in the main). However in the spin-independent amplitudes 
two sets of LEC's enter. One set is the pair of LEC's which contribute to a and /3; these have 
been fixed by requiring the full polarisabilities (with third and fourth order contributions) 
to take the experimental values. The other set are the parameters Cj from £( 2 \ which enter 
in the tadpole graphs. These are rather poorly known; furthermore they have been deter- 
mined from pion-nucleon scattering at second, third and fourth order with quite different 
results (and indeed from different data sets with different results). Clearly the second-order 
determinations are the most appropriate to use in these calculations as vertex dressing is 
absent at this order. We have used c\ = —0.81 GeV, c 2 = 2.5 GeV and c 3 = —3.8 GeV |2~9~ 
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Fig. 5: Effect of altering the spin polarisabilities by hand. Solid line: 4th order HBCPT 
(as in previous figures); long dashes, short dashes and dots: polarisabilities taken from the 
DR analyses of refs JI3| , |I5| , |n| respectively. 



Sensitivity to these choices is explored later. 

All our amplitudes are calculated in the Breit frame. At third order, the result is the 
same in any frame (so long as the energies remain low); the amplitudes could be in the Breit 
frame, or the centre-of-mass frame, or even the lab frame (since uj — uj' is of higher chiral 
order). At fourth order however the result will be different in different frames. Interestingly 
there is still no dependence on uj — uj', but there will be expressions involving the average 
nucleon three-momentum p. In the Breit frame p vanishes, so the amplitudes take on by 
far their simplest form. In the centre-of-mass frame p can be written in terms of the photon 
momenta, and the structures which enter are unchanged, but the amplitudes will have terms 
with the "wrong" symmetry in uj. (These of course can be obtained from a boost of lower 
order amplitudes, and are not independent.) The differential cross section in any frame is a 
kinematic prefactor $ 2 times the invariant |T| 2 , which in terms of the Breit-frame amplitudes 
is given by 

|T| 2 = \A\(l + cos 2 B) + |A 2 (3 - cos 2 9) 



+uj 2 sin 2 9 [4v4 3 v4 6 + (A 3 A 4 + 2A 3 A 5 - A^ 2 ) cos#] 
+uj 4 sin 2 9 \\A\ sin 2 9 + \A\{\ + cos 2 9) + A 2 (l + 2 cos 2 
+2A 6 (A 4 + 3A 5 ) cos 9 + 2A A A 5 cos 2 9] 



3AI 



(3) 



Above the threshold one should replace A\ by |A| 2 and AiAj by §l(A*Aj). 

Experimental data, though taken in the lab frame, are sometimes presented as centre- 
of-mass frame differential cross sections, so both prefactors are needed: 



1 J 

$iab = - ; 

47T UJ 



1 m N 

47T y/s' 



(4) 



There is a further complication which stems from the heavy baryon reduction. Physically, 
there is a threshold at a lab energy of uj^ = M n (l + M 7T /2m N ). However the HBCPT 
amplitudes at third order have cusps at uj = M w , and at fourth order they are actually 
singular at this point. The resolution was given by Bernard et al. when discussing forward 
scattering M: if calculated to all orders, the sum of all the singular terms will be precisely 
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what is required to shift the threshold to the physical value. Thus at finite order, to have 
the threshold at the right point and maintain finite amplitudes, one changes variables from 
the Breit frame energy uo to ( = cu/c^th, then Taylor expands and discards pieces of higher 
order. At fourth order, therefore, if a naive loop amplitude is f i3) (u/M n ) + fW(u/M n ), the 
shifted amplitude is 

The subscript "LO" indicates the dropping of terms of 0{l/m 2 N ). The shifted amplitude is 
finite at ( = 1. 

Since |T| 2 is an invariant, it should make no difference in which frame it is calculated, 
and we use the Breit frame as the amplitudes are simplest there. Of course in HBCPT 
the invariance is only respected up to the order in l/m N to which we are working, and so 
the results obtained by calculating the amplitudes in different frames would not in fact be 
identical. 

Although we have calculated the amplitudes consistently to fourth order, those ampli- 
tudes substituted in Eq. |3] will not give a differential cross section which is of consistent 
order. Since A\ starts at second order with the Thomson term, knowledge of the amplitude 
to fourth order gives the cross section to sixth order. The fourth order parts of all the 
other Ai are not actually needed, and the seventh- and eighth-order parts of \Ai\ 2 should 
be subtracted. We do not do this. The various orders of contributions are actually of very 
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Fig. 7: Effect of varying the pion-scattering LEC's q. Values are taken from Fettes and 
Meifiner 291. Solid and dashed lines are the second order fits based on KA85 and SP98; 



dotted is the fourth order fit to SP98 from Table 1 of that work. 



similar magnitude, and if we try to remove higher order pieces we may not even be left 
with a positive cross section. Thus all the cross sections we present are calculated with 
fourth-order amplitudes, and our kinematic prefactors and change of variable from u to ( 
and from one frame to another are done exactly. 



III. RESULTS 

In Figs. 2 and 3 we show differential cross sections at various lab and centre-of-mass 
angles for which data exist below 200 MeV. In each graph we show the Born (pion and 
nucleon pole) contribution at fourth order, and the HBCPT results with the amplitudes 
calculated to third and fourth order. The third as well the fourth order results have had 
the variable change u/M^ — > uj/uj^ as explained above; otherwise the threshold would be 
in the wrong place. This, and the fact that we take the basic third-order amplitude to be 
in the Breit frame, rather than the centre-of-mass frame, accounts for the slight differences 
between our third-order results and those of Babusci et al. ||20]| . In Fig. 4 we show three 
third order curves: amplitudes calculated in the Breit frame and shifted, in the Breit frame 



and unshifted, and in the CM frame and unshifted as in ref. ||20|| . The differences between 
these are higher order. 

For scattering angles up to 90° the fit at either third or fourth order is very good except 
at the highest energies, and much better than the Born terms alone. In the cusp region 
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the fourth order is arguably better than the third order at all angles. The good fit to the 
data above threshold is a feature which has not been remarked before. However for larger 
angles both curves start to fall below the data as the threshold is approached, and fall 
dramatically short above. Of course one does not expect agreement to continue indefinitely, 
as in experiments the A resonance is clearly seen, and this will not be reproduced by HBCPT. 
However it is not obvious why backward angles should be so much worse than forward angles. 
As the data in that region is only from one experiment [EJ one might question its accuracy, 
but in fact it fits well with dispersion analyses which incorporate the more numerous higher 
energy data. 

At fifth order, the new diagrams which enter fall into several groups. First, there are 
two-loop graphs, which one might well expect to be small. Then, there are one-loop graphs 
with two insertions from C^ 2 \ Some of these give higher l/m N corrections and these again 
ought to be small (some indeed have already been incorporated in the change of variables 
which corrects the threshold) but some, such as two insertions of the anomalous isovector 
magnetic moment, may be numerically significant. Then there are one- loop graphs with 
one insertion from C^ 3 \ and seagulls contributing directly to the polarisabilities. It is in 
these that the effects of the A will show up in the spin-dependent amplitudes, and these 
contributions are expected to be anomalously large. We have not calculated fifth-order loop 
graphs. However the seagulls may be incorporated effortlessly, since the free LEC's will just 
be fit to give the "experimental" values of the 7$. In Fig. 5 we show the effect of setting the 
7i to the values found in three fixed-t dispersion-relation studies Hi3Hl5|j. Below threshold 
the agreement is slightly improved, but above the enhancement at backward angles is not 
sufficient to give agreement. 

Finally, we show the effects of varying the input parameters within a purely fourth-order 
calculation. These are the places where the effect of the A is already present; for instance 
the tadpole graphs of Fig. ln-r could arise from integrating out an intermediate A. These 
graphs vanish for forward scattering but contribute significantly to the spin-independent 
amplitudes at backward angles. First we show the results of varying a — (3 while keeping 
the sum fixed. This shows up only at backward angles; to reproduce the Mainz data at 180° 
we need to make a dramatic shift, to a = 8 and (3 = 6 x 10~ 4 fm 3 . (Rather similar numbers 
were found by GrieBhammer and Rupak in a study of Compton scattering from the deuteron 
in the framework of effective field theory without dynamical pions [0.) However even these 
extreme values do not cure the problem above the threshold. Then we show the effect of 
using a very different parameter set for the q that enter the tadpole graphs, while keeping the 
polarisabilities fixed. (These graphs do of course contribute to the polarisabilities, as given 
ref. 0. However the fourth order LEC's are adjusted to give the measured polarisabilities 
whatever the other contributions might be.) These parameters are quite uncertain and differ 
substantially in different fits; we show two second-order fits and one fourth-order pS], Again, 
the only effect is at backward angles, and none of the parameter sets gives agreement above 
the threshold. 

Thus the conclusion of this fourth-order study is similar to that of the third-order study 
of Babusci et al. f2(|. For angles up to 90° , HBCPT without an explicit A does an excellent 
job of fitting the experimental differential Compton scattering cross section up to about 
170 MeV. However it fails at larger angles for photon energies above about 120 MeV. Very 
large fifth-order contributions would be needed to approach the experimental data above 
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the threshold at these angles. 
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APPENDIX A: FOURTH-ORDER LOOP AMPLITUDE IN THE BREIT FRAME 



The full amplitudes in the Breit frame for the diagrams of Fig. 1 are as follows. The 
notation tj is used for the tensor structures which multiply the amplitudes A4 of Eq. [1]; for 
example t± = e' ■ e. 



T 



9 2 e 2 

^m N f 2 



(t 1 + t 3 ) (d-i) 



dJ 2 [u, m? 



m f 1 r 

T b = ~o f2 / dX 

2«W; Jo 



(t! + t 3 )(d+l) 



dJ 2 [xu, m 2 ] 



2 1 -,dJo[u!,m 2 



duo 



+crossed 



dxu 



/ 9/ 1 \ ,00 in/ n\ dJnlxu, m 2 ] 1 
■[xuj 2 (t 2 + it 5 ) + {x w 2 + lxt)(tx + t 3 )J — + crossed 



T, 



g 2 e 2 



2m N fl 



■r 3 (ti + t 3 )(ujj [uj,m 2 } + A^-fm 2 ]) + crossed 



9 2 e 2 
2m N f 2 



r 3 (ti + t 3 ) J dx(u Jq[xu, m 2 ] + A,j[m 2 ]j + crossed 

g 2 e 2 1 
T e = — (1 - r 3 )(ti + t 3 )-(j 2 [uj,m 2 } - J 2 [0,m 2 ]) + crossed 



g 2 e 2 n 



4m N f 2 



J dx (fj, v — /i s r 3 ) (((1 — 2x)t 3 + ti) cos 9 — t 2 — (1 — 2x)t 4^ uj J $[xoj. 



m 



-(1 — r 3 ) (2(d + l)tiu; 1 J 2 [xu, m 2 } + (2x(l — x)ti — |(1 — 2x)t q)uj J o[xu, m 2 ] + crossed 



9 2 e 2 



g 2 e 2 n 



([/jl v + /i s r 3 )((ti + t 3 ) cos6> - t 2 + t 4 - t 5 ) + |(1 + r 3 )t 6 ) u dxJ [. 

dJ 2 [u, m 2 — xyt] 
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+ crossed 



. dy dx 
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duj 



+ (ti(2(d + 3)V(x, y)-t)+ uj 2 t 2 ((d + 3)((d + 5)xy - 2x - 2y) + 2) + 

((d + 3)x - lV 2 (t 6 - t 5 ) + (d + 3)(l-x-y V 2 t 4 )— m ' " ^ 
+ (-t 1 V(x, y) 2 - (art 6 - arte + (l-x-y)t 4 ) y)w 2 



<9cD 
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crossed 
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(d 2 - l)A^[m 2 - x(l - x)t\ + 2(1 - 2(d + l)x(l - x))tA' n [m 2 - x(l - x)t\ 
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+ ((2 Cl M 2 - (c 2 + c 3 )^ 2 + V(x, y)c 3 )ti + 

(xyc 2 + ((d + 4)xy - x - y)c 3 )cj 2 t 2 ) A^m 2 - xyt] 
+xyuj 2 t 2 (-A Cl Ml + 2w 2 (5 2 + c 3 ) - 2F(x, y)c 3 ) A^[m 2 - xyt]" 
2e 2 



f 2 

J 7T 



dx 



(dc 3 + c 2 )A n [m — x(l — x)t] 



■(4 Cl M 2 + 2x(l - x)tc 3 )A;[m 2 - x(l - x)t] 



2„2 



We 
8m N f 2 



;i + r 3 )(d-l)t 1 A 7r [m 2 



(Al) 



where the q are low energy constants from £^1, and c 2 = c 2 — g 2 /8m N . The integrals 
J [o;,m 2 ], J 2 [cj,m 2 ] and A T [m 2 ] have their usual meanings, prime denotes differentiation 
with respect to m 2 , u = (1 — x — y)u and 



V(x, y)=u 2 + \t{\ -x-y + 2xy). 



(A2) 



We have also introduced an extra tensor structure, tr — cr • (q' x q)e' • qe • q'. This is not 
independent; t 7 = sin 2 #t 3 + cos#t 5 — t 6 , but as it arises naturally in the calculations — and 
enters at too high an order in uj to affect the polarisabilities — it is a useful notation. 

The notation "+ crossed" means that to every term, another is added with e <-> e', 
q <-> — q' and uj <-> —a;. Since the tj are all either symmetric or antisymmetric under this 
transformation, the net effect is to add a term with uj <-> — c<j to the coefficients of ti and t 2 , 
and subtract such a term from the coefficients of t 3 . . . t 7 . 
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The spin- dependent terms do not all look identical to the expressions given in ref. 
but they are equivalent. 
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